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Abstract 

General expressions are given for the generation of Primordial Black 
Holes (PBH) in a universe with a presently accelerated expansion due 
to a(n effective) cosmological constant. We give expressions both for a 
powerlaw scalefree primordial spectrum and for spectra which are not 
of that type. Specializing to the case of a pure cosmological constant 
A and assnming flatness, we show that a comological constant with 
Ga,o = 0.7 will decrease the mass variance at the PBH formation time 
by abont 15% compared with a critical density nniverse. 
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1 Introduction 


The generation of a spectrum of primordial fluctuations in the very Early Uni¬ 
verse is the crucial ingredient of all inflationary scenarios. These fluctuations 
can explain the generation of all (classical) inhomogeneities that can be seen in 
our universe, from the Cosmic Microwave Background (CMB) anisotropies to the 
Large Scale Structures (LSS) in the form of galaxies and clusters of galaxies. The 
inflationary paradigm therefore reconciles Big Bang cosmology with the appear¬ 
ance of an inhomogeneous universe |]l|. In addition, each inflationary scenario 
makes accurate predictions allowing for the observations of ever increasing vari¬ 
ety and quality to discriminate between the various model candidates. One such 
prediction is the possible formation of Primordial Black Holes (PBH). Indeed, it 
was realized already some time ago that a spectrum of primordial fluctuations 
would lead to the production of PBH |Q. For this generation mechanism to be 
efficient, one typically needs a “blue” spectrum [^. In this way, one can hope 
that the density contrast averaged over the Hubble radius is sufficiently large 
that the resulting PBH production is not unsignificant and can be used as a 
powerful constraint on the spectrum of inflationary primordial fluctuations and 
the underlying high-energy physics model [|, |^. The production of PBH takes 
place on scales much smaller than those probed by the CMB anisotropy and LSS 
formation. In this sense, it is analogous if less spectacular, to the generation of 
a primordial gravitational wave background in inflationary models. Of course in 
the latter case, its discovery would be a remarkable prediction of inflation while 
the existence of PBH is a confirmation of the existence of the primordial fluctua¬ 
tions spectrum itself, irrespective of the way it was generated. In a recent paper 
1^, it was shown that the mass variance unitk) at the PBH formation time tfc, 
was significantly overestimated due essentially to an incorrect relationship be- 
tween anitk) and Snitk) or equivalently ^(k,tk) which quantifies the power 
at Hubble radius scale (these quantities will be carefully introduced below). This 
must be corrected if one is to make accurate predictions. In 0, only a critical 
density universe was considered and expressions were derived for scalefree pow- 
erlaw spectra. Recent Supernovae observations strongly suggest that we live in a 
presently accelerating universe with ~ 0.3, r2A,o ~ 0.7, the inclusion of a(n 
effective) cosmological constant seems further to make all observations converge 
into a consistent picture. Hence, it is important to derive general expressions in 
our formalism valid for a flat universe with an effective cosmological constant like 
in quintessence models. In the particular important case of a pure cosmological 
constant A, also considered in 0, the dependence on A can be quantified accu¬ 
rately and the relative decrease we find is in agreement with part of the analysis 
done in 0]. Finally we also generalize our results to primordial perturbations 
spectra which are not scalefree, an interesting possibility to consider in view of 
the wide range of scales probed by PBH formation. We first review the formalism 
describing PBH formation. 
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2 PBH formation 


We assume for simplicity that a PBH is formed when the density contrast av¬ 
eraged over a volume of the (linear) size of the Hubble radius satishes 5min < 
^ < ^max, and further that the PBH mass, Mpbh, is of the order of the “horizon 
mass” Mh, the mass contained inside the Hubble volume. Relyimg on semiana- 
lytic considerations it is common to take 6min = Smax = 1 but recent numerical 
simulations suggest rather 6min ~ 0.7 [@| and show that Mpbh can span a certain 
range, around Mp though, at a given formation time. More accurately, when 
some scale dehned by the wavenumber k reenters the Hubble radius after infla¬ 
tion at some time with k = {aH)\t^, it can lead to the production of PBH 
with MpBH ~ Mpitk)- Obviously, there is a one-to-one correspondence between 
and k. 

For Gaussian primordial fluctuations, the probability density Pr{S), where 6 
is the density contrast averaged over a sphere of radius R, is given by 


Pr{S) 


1 _ 

_ g 2cr^ (K) 

\/^ ciR) 


( 1 ) 


Here, the dispersion (mass variance) cr‘^{R) 
top-hat window function. 



is computed using a 
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a\R) = ^J dkk^ W^HikR) P{k) , 


( 2 ) 


where P{k) is the power spectrum (we assume isotropy of the ensemble). From 
a point of view of principles, the averages are quantum averages; however, an 
effective quantum-to-classical transition is achieved during inflation . For PBHs 
produced by inflationary perturbations, this quantum-to-classical transition is 
guaranteed for all masses of interest to us (see |T0[|). 

The expression WpHikP) stands for the Fourier transform of the top-hat 
window function divided by the probed volume Vw = ^tiP^, 


WpHikR) = 


{kpy 


(sin kP — kP cos kP) . 


( 3 ) 


Hence the probability (3{Mh) that a region of comoving size R = 
an averaged density contrast at horizon crossing tk in the range 6min < ^ < 5max, 
is given by 


P{Mh) 






V^27r 6min 


;2 


mzn 




( 4 ) 


where crjj{tk) = a‘^{R)\^^, and the last approximation is valid for 6min S> crpitk), 
and i^Smax Antin') O'Hitk')- 
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Important conclusions can be drawn from (^). Let us consider first the value 
of (3{Mh) today. Today we have cT|^(to) — 10“®, so clearly the probability of 
forming a black hole today is extraordinarily small. This probability can increase 
in the primordial universe if the power is increased when we go backwards in 
time, but the probability will remain very small, j3{Mii) -C 1, at all times due to 
the magnitude of > 1. 


3 Mass variance in the presence of A 

When the universe contains a cosmological constant A, this must be taken into 
account for a correct accurate calculation of the mass variance at early times. In 
this section, we will extend the formulas derived in p. As stressed there, one 
should distinguish the behaviour of the quantity crjjitk), which is ultimately the 
quantity of interest, from the quantity k^(f)‘^{k,tk) or 6jj{k,tk) with 

Snik^tk) = ^ P{k,tk) = , ( 5 ) 

where tk is the PBH formation time of interest, deep in the radiation dominated 
stage. However, when dealing with a flat universe with < 1, we have today 

/p3 2 

to) ^ ^ P(fco, to) = n-J^o kl ^\ko, to) , (6) 

where /o stands for any quantity / evaluated today (at time to), ^m,o = 
is the present energy density of dust-like matter relative to the critical density 
and r 2 A,o = 3 ^ = 1 — ^m,o- We hrst relate the quantities appearing in (H,H) at 
the formation time tk and at the present time to for arbitrary evolution of the 
universe after radiation domination and for a scalefree powerlaw spectrum. 

General expressions with powerlaw spectrum: assuming a scalefree power- 
law primordial spectrum of the type k^^’^{k) = A{t) k'^~^ on super Hubble radius 
(“superhorizon”) scales, we then have 

k^ ^^{k,tk) = (^ (^l-^j^adt'^ kl ^^{ko,to) , (A 

and analogously 

SUkAk)=(^ (^l-^J^adt'^ ^nikoAo) ■ ( 8 ) 

The lower limit of integration in can be safely taken to be zero. In (|^,|), 
we have used a radiation dominated stage followed by some arbitrary evolution 
of the scale factor. In earlier work, we considered a radiation dominated stage 
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followed by a matter dominated stage which constitutes a special case of 
An effective cosmological constant as in quintessence models would also be a 
particular case of (0J|). However, in contrast to a pure cosmological constant 
A, the time evolution of the scale factor a{t) is model-dependent and cannot be 
given in full generality. The quantity fcg <h^(fco,to)) or equivalently Sjj{kQ,to), at 
the present Hubble radius scale can be derived using the large angular scale CMB 
anisotropy data. It is that quantity that comes from observations which hxes the 
overall amplitude of the fluctuations spectrum. The COBE data show that it is 
of the following order of magnitude 
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q 4^(4,*o) =0.86 X 10 ® /lj({nj}) , 


( 9 ) 


where A^{{ni}) parametrizes the amplitude variations and is chosen such that 
Alin = 1, = 0.3) ~ 1 Alin = 1, = 1) - 1-94 . (10) 


The exact amplitude depends on the cosmological parameters {rij}, referring to 
the background as well as to the inflationary perturbations, and this model de¬ 
pendence is encoded in the quantity of order unity Ao({nj}). Eq. (|l0|) assumes 
a powerlaw spectrum with spectral index n at least on large scales. For hxed 


n 7 ^ 1, while the absolute values in (pTl) are modihed, the ratio between them is 


unaltered []^. Note that for a quintessence model, Ali{ni\) is model dependent. 

Finally, we must relate all our results to the quantity of interest for the com¬ 
putation of the PBH abundance, the mass variance crjjitk) on the Hubble radius 
scale at horizon crossing time As stressed in 0, one has (with k = iaH)\t^) 


^ (fc) S^ikAk) 


( 11 ) 


It is crucial to distinguish both quantities auitk) and (5_H-(fc, 4). As seen from (Q). 
the quantity crjjitk) which depends on the averaging procedure through eq.(^ is 
correctly related to the (non-averaged) quantity Sjjik,tk) in a non-trivial way. 
The quantity Snikytk) can be reconstructed at the time tk from its present value 
Snikoyto) using (H). But this is not the case for the quantity anitk) because the 
deformation of the power spectrum is different at the time tk and today. In other 
words, 

Tik\to)^Tik'Ak) , (12) 

where the transfer function T(fc,t) is dehned through 

Pik, t) = Pik, U) T^ik, t) , T(fc ^ 0, t) ^ 1 . (13) 

Here, ti is some initial time when all scales are outside the Hubble radius, k <C aif, 
we can take for example ti = t^, the end of inflation. For a powerlaw scalefree 
spectrum we have 0 

fee 

a^ik)= I T‘^ikx,tk) W^fjix) dx , (14) 

Jo 
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where ke corresponds to the shortest fluctuations wavelengths with the size of 
the Hubble radius at the end of inflation. The transfer function T{k',t) in the 
integrand of (JT^) must be taken at the time tk, not today. The accurate value of 
a{k) requires numerical calculations but estimates made in |p show a significant 
overestimation of the mass variance (Tnitk) when it is not computed correctly 
using the right quantity a{k). One then gets that (10/9)^a^(fc) -C 25 (the value 
usually taken in the literature for a critical density universe) for all the mass range 
of PBH produced in the radiation era. If one is willing to use PBH formation as 
a precision tool in cosmology, it is important to check in how far the presence of 
a cosmological constant with r2A,o = 0.7 brings further modifications. 

An important conclusion can be immediately drawn by inspection of the in¬ 
tegrand in ([T^) without accurate knowledge of the transfer function at time tk- 
Indeed, as we are interested in times tk -C 4? and in universes where Ha domi¬ 
nation occurs late, it is clear that at the time tk, neither the long-wave nor the 
short-wave fluctuation modes are affected in any way by the presence of (an ef¬ 
fective) A. As can be seen from (^), this implies that the transfer function at 
time tk does not depend on A and the same must apply therefore to the quantity 
a{k). We conclude that the influence of a(n effective) cosmological constant on 
the probability P{Mh) comes solely from its influence on the quantity 5‘jj{k,tk), 
or k^ <f)^(fc, tfc). It is this influence that we will quantify in the next subsec¬ 
tion. We now consider a powerlaw spectrum and specialize to a universe with a 
cosmological constant A. 


Powerlaw spectrum with A: in order to account for the presence of a cos¬ 
mological constant A, we must replace the evolution of the scale factor a{t) after 
the radiation dominated stage. The scale factor for this stage of the universe 
evolution is very well approximated by 0 


a{t) = ai sinh^'^^(/3t) , 


(15) 


where |/3 = -^/A/S = TToa/HaO) =const. The evolution (|T^) smoothly inter- 
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polates between a pure (flat) dust-like matter dominated stage, with a{t) oc fa, 

for -C 1 which is of course the case right after teq, and a A dominated universe 

in the asymptotic future. In particular for scales for which z^q > z{tk) 3> 1, 
2 

a{tk) oc tl- It is this evolution ([T5| ) which must be used in (|^,|[). It is physi¬ 
cally appealing to express the results in terms of the quantity Mnitk)- Then the 
following result is obtained 


k^ ^^{k,tk) = 


X 



k 1 

^eq 

n—1 

Mniteq) 

ko _ 


_ Muitk) _ 




( 16 ) 


The evolution ( 0 ) must now be substituted in (pT]). Actually it is slightly more 
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( 17 ) 


accurate to write the quantity ^ as 

k 1 

— cx MnitkY^ , 
ko 


tk -c t 


eq 


where the proportionality constant is independent oi (and hlA,o)- We will take 
now the following quantities: = 9.81 x 10“® relative 

energy density of relativistic matter today, peff is the present effective number of 
relativistic degrees of freedom, = -^km/s/Mpc, T.^_o is the present tempera¬ 
ture of the CMB. Using these quantities, eq. ([T6|) hnally becomes for IIa.o = 0.7 


= 1.75x10" 


0.219 ,2 


/2(flA,G 




X 


9.75 X 10^® k 


65 


O 

n—1 

g 

/ 2.726°K 


Mnitk) _ 


n —1 
2 


.( 18 ) 


where Muitk) is expressed in grams. In deriving (pISD we have used 


Ha 


f*to 


J(r2A,o) = 1-/ a(f) df ~ 1 


no Jo 


3y/flA,0 


r 

sinh X 

lo 

sinh d 


, 1 1 -|- ■\/flA ,0 

a = - In ■ 


2 1 — ■\/hlA ,0 

which gives the following numerical result substituted in 

/(0.7) = 0.468 . 


dx (19) 

( 20 ) 

( 21 ) 


Hence, /“^(0.7) = 4.57 gives an increase of 64% compared to the value 25/9 = 
2.78 obtained for a critical density flat universe, flm.o = 1- On the other hand 
normalization to the CMB fluctuations decreases Hq by a little bit more than 48%, 
as seen from ( 0 )- Therefore both effects combined lead to a decrease of about 
15%. As for a{k), we have seen above that it is insensitive to the presence of a 
cosmological constant. From (^), the correspondence between the (approximate) 
PBH mass M^itk) and k is independent of for tk "C tgg- Therefore we 
conclude that, in a flat universe, a cosmological constant with Ha.o = 0.7 will 
decrease the mass variance (Tnitk) as follows 

^1(4)Ioa, 0=0.7 - 0.85 a|(4)|o^,o=i (22) 

As a result, the signihcantly lower value estimated in [^] is further reduced by 
about 15%, further diminishing the probability for PBH formation. This is in 
good agreement with the analysis performed in [0. Still, it is instructive that 
these authors have incorrectly assumed that anitk) = C duitk) where the con¬ 
stant C is the same at all times tk up to the present time to- Nevertheless, their 
eq.(4.30) would give essentially the same relative decrease of anitk) due to the 
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presence of a cosmological constant as in (^). This is because the correct factor 
a{k) does not depend on A. 

We hnally note that equation (p!8[) makes use of the observed amplitude today 
on the present Hubble radius scale and, as far as perturbations are concerned, 
combines it with an assumed (powerlaw) behaviour towards small scales on a very 
broad range of scales. Other behaviours are certainly possible however. This is 
reminiscent of the primordial gravitational wave background generated during 
inflation extending up to frequencies as high as lO^^Hz. There too, one can 
imagine a behaviour towards large frequencies departing from a simple scalefree 
law (see e.g. []A| for such a model with a jump in the tensorial spectral index ht). 


Hence if the assumption of a scalefree spectrum does not hold, a more general 
expression will be needed. For this reason, we now generalize our results also to 
these cases. 


Spectrum with a characteristic scale: A further important generalization 
concerns the primordial fluctuations spectrum itself. Indeed, the equations writ¬ 
ten in the previous subsection assume a scalefree spectrum. However, this needs 
not be the case especially in view of the large range of scales that are probed by 
PBH formation. 

Let us therefore dehne in full generality 

=(l--[ a dA F{k) = ^ F{k) , (23) 

V a io Jt=t, 9 

where F{k) can be any complicated function of k. The expression ([^) represents 
the primordial spectrum on “super-Hubble radius” (superhorizon) scales. For 
example, the spectrum of double inflation considered in |^, is of this general 
type. The corresponding generalization of (|^), or (^, leads to a more complicated 
equation, viz. 


k^ <F^(fc,4) = 


1 - 


H 


-2 


a dt' ) /cq <F^(A;o, to) 


t=to 


X 


F(aiM^"(4)) 

F{ko) 


(24) 


where ai=constant. In case the function F{k) is actually of the form F"(^), 
where kg dehnes the characteristic scale - an example of such a spectrum was 
found in and considered in - eq.(|^) can be recast into a slightly simpler 
form 


k^ ^^{k,tk) = 


H w 

— / a dt 


-2 


t=to 


kl <F^(fco,to) 


F 


X 


Mnitk) I 


F{ko/ks 


(25) 
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with Ms = MnitkJ and tk^ < teq- Specializing to the particular case = 0.7 


just requires the substitution, like in 


of the corresponding numbers into 


Finally, we come to the calculation of unitk) itself. We now have the corre¬ 
sponding generalization to primordial spectra of arbitrary shape 


a 


H 


(4) = 


8 


Slvr^ 


F{kx) T'^{kx,tk) W^u{x) dx 


In particular, the general expression for a^{k) is given by 


a 


\k) = 


F{kx) 

'm 


x^ T‘^{kx,tk) W^j^{x) dx . 


(26) 


(27) 


Again, for the case of interest to us, tk -C teq, the transfer function at the time 
tk and therefore also a{k) are independent of A. We conclude that the same 
decrease found in (P^ ) will apply here too. Note that (^) extends the result 
(0) derived in [Q for a powerlaw scalefree spectrum which just corresponds to 
F{k) oc k'^~^. A characteristic scale in the primordial spectrum is an interesting 
possibility with respect to PBH formation in view of the large range of scales 
involved, much larger than CMB anisotropy or LSS formation. We have already 
considered some simple toy models in |P and interesting results were obtained. 
The detailed numerical study of more sophisticated spectra and their possible 
relevance to observations is under progress 
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